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We theoretically study the electronic band structure of isolated unpaired and paired dangling 
bonds (DB), DB wires and DB clusters on H:Si(100)-(2x 1) surface using Extended Hiickel Theory 
(EHT) and report their effect on the Si band gap. An isolated unpaired DB introduces a near-midgap 
state, whereas a paired DB leads to n and 7r* states, similar to those introduced by an unpassivated 
asymmetric dimer (AD) Si(100)-(2xl) surface. Such induced states have very small dispersion due 
to their isolation from the other states, which reside in conduction and valence band. On the other 
hand, the surface state induced due to an unpaired DB wire in the direction along the dimer row 
(referred to as [110]), has large dispersion due to the strong coupling between the adjacent DBs, 
being 3.84A apart. However, in the direction perpendicular to the dimer row (referred to as [110]), 
due to the reduced coupling between the DBs being 7. 68 A apart, the dispersion in the surface state 
is similar to that of an isolated unpaired DB. Apart from this, a paired DB wire in [110] direction 
introduces tt and 7r* states similar to those of an AD surface and a paired DB wire in [110] direction 
exhibits surface states similar to those of an isolated paired DB, as expected. Besides this, we report 
the electronic structure of different DB clusters, which exhibit states inside the band gap that can 
be interpreted as superpositions of states due to unpaired and paired DBs. 

PACS numbers: 73.20.-r, 73. 20. At, 74.78. Na 



I. INTRODUCTION 

Due to their importance in Si technology, states in- 
duced by the dangling bonds (DBs) on Si(100) surface, 
known as P5 centers, have been a topic of study for 
decade^^^^ii^^^^. These are important 
due to their role in introducing states inside the band 
gapi£ i n i 12 i 13 i 16 , Fermi level pinning^ 7 -, charging issue - 
consider variation of threshold voltage due to charged 
DBs, reliability issues - consider negative bias tempera- 
ture instability 18 , enhancement of 1/f noise due to trap 
states resulting in poor DC characteristics^. With ad- 
vances in nanoscale science, detailed understanding of 
these surface states is needed to exploit their role in pos- 
sible technological applications, e.g. use of DB as a tem- 
plate for the growth of molecular nanostructures 2 ^. This 
realization has been the motivation behind research like 
desorption of H from Si surface forming a D B 21 i 22 i 23 , elec- 
tronic and electrostatic effect of DB on transport through 
a styrene wir o 16 i 24 , DB dynamics on Si surface 25 , Jahn- 
Teller distortion in DB wires at low temperature^, mag- 
netism in DB structure s 27 ! 28 ! 29 , passivation of DBs with 
H and D 30 , etc. 

Watanabe et al 11 predicted that an unpaired DB (Si 
dimer with one hydrogenated atom) wire on H:Si(100)- 
(2x1) surface in the direction along the dimer row would 
introduce states inside the band gap. This prediction was 
later experimentally observed, however the mechanism 
is still not clear—. Recently, Cakmak et aJ^l theoreti- 
cally reproduced their results and found similar trend. 
All these calculations were based on Density Functional 
Theory (DFT) where band gap is underestimated 3 - 2 .. It 
is thus hard to quantitatively benchmark the results. In 
this paper, we use Extended Hiickel Theory (EHT) for 



electronic band structure calculations using transferable 
parameters developed by Cerda et a&. Notably, EHT 
gives correct Si bulk band gap and other features as 
reported by Kienle et aJ2i. It has been applied to Si 
surfaces obtaining satisfactory result s 34 ! 35 . The state of 
the art in electronic structure calculations for Si is GW 
approximatio n 36 ! 37 . However, computational complex- 
ity prohibits its use in transport calculations for large 
systems like the ones discussed in this paper where the 
Si atoms per unit cell may approach 384. EHT pro- 
vides computationally accessible, yet accurate approach. 
Apart from this, to the best of our knowledge there has 
not yet been a systematic study of isolated DBs, DB 
wires and DB clusters on H:Si(100)-(2xl) surface using 
a method that quantitatively addresses the surface states. 

In a previous study-i^, we report that a DB electroni- 
cally affects Si atoms up to approximately 10 A away from 
it, by introducing a near-midgap state in the local density 
of states of the neighboring Si atoms. Therefore, to simu- 
late DBs in isolated, wired and clustered configurations, 
we use a large enough unit cell in the directions parallel 
to and perpendicular to the dimer row, conveniently re- 
ferred to as [110] and [110] directions respectively 3 - 2 - (as 
shown in Fig. 1(e)). Few large unit cells are shown in 
Fig. 2 to demonstrate the approach. This approach has 
been used previously in Ref j n i 16 ! 31 . Table I summarizes 
different surfaces being investigated. We use electronic 
band structure principles in [110] and [110] directions. 
However in [100] direction, i.e. direction normal to sur- 
face, we simulate a semi-infinite structure by using six- 
teen layers of Si atoms as shown in Fig. l(a-d). The im- 
portant aspect of this approximation is the use of a large 
enough unit cell in [100] direction, which minimizes the 
quantum confinement effects and the correct Si band gap 
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FIG. 1: (color online) Ball and stick model of sixteen layer 
Si(100)-(2x 1) unit cell used to construct different surfaces. 
Top four layers, represented by J, are relaxed due to surface 
reconstruction. Bottom twelve, represented by f, are bulk 
layers, (a) Unpassivated asymmetric dimer (AD) unit cell re- 
ferred to as paired dangling bond (DB) configuration, (b) Un- 
passivated symmetric dimer (SD) unit cell, (c) Hydrogenated 
SD unit cell, (d) Hydrogenated unit cell with one DB referred 
to as unpaired DB. Arrow shows [100] direction, perpendicu- 
lar to the surface, (e) H:Si(100)-(2x 1) surface with direction 
along and perpendicular to dimer row labeled as [110] and 
[110] respectively. In part(a-e), back surface is hydrogenated 
to eliminate any DB induced states, (f) A SiqHia cluster with 
H-Si bond length of 1.47A and H-Si dimer angle of 111°. 



for H:Si(100)-(2xl) is obtained (as shown in Fig. 3(a)). 
Fig. 3(a) also shows the tt and 7r* states for unpassivated 
asymmetric dimer (AD) and unpassivated SD surface, re- 
ferred to as Si(100)-(2xl)-AD surface and Si(100)-(2x 1)- 
SD surface respectively. Fig. 3(b) shows transmission 
through the above-mentioned three surfaces. H:Si(100)- 
(2x1) surface has a clean band gap, Si(100)-(2x 1)-AD 
surface has 0.62eV of band gap and Si(100)-(2xl)-SD 
surface has no band gap. 

We present calculations of isolated unpaired DB (Si 
dimer with one atom hydrogenated) and paired DB (un- 
passivated Si dimer) in Fig. 4(a) and (b) respectively. 
An unpaired DB introduces a near-midgap state as we 
previously report^, whereas a paired DB introduces a 
7r state close to the valence band maximum (E v ) and a 
7r* state close to the conduction band minimum (E c ). 
These states due to a paired DB are similar to the tt 
and 7r* states introduced by Si(100)-(2x 1)-AD surface 
as shown in Fig. 3(a), however dispersion is smaller than 
that of AD surface states. We provide details of band- 
widths of different states in Section IV. Furthermore, we 
model DB wires consisting of unpaired and paired DBs 



TABLE I: Details of the unit cells used to simulate different 
surfaces. Sixteen layer atomic structures shown in Fig. 1 are 
used to construct bigger unit cells. Thus, a (4x6) unit cell has 
dimensions equal to 15. 36 A and 23.04A in [110] (perpendic- 
ular to dimer row) and [110] (along dimer row) directions re- 
spectively. The total number of Si atoms (excluding H atoms 
on top and back surface) per unit cell are given. 



Unit cell surface 


Unit cell a 


Atoms 6 


H:Si(100)-(2xl) 


2x1 


32 


Si(100)-(2xl)-AD 


2x1 


32 


Si(100)-(2xl)-SD 


2x1 


32 


Unpaired DB 


4x4 


256 


Unpaired DB wire 


4x1 


64 


Unpaired DB wire d 


2x4 


128 


Paired DB 


4x4 


256 


Paired DB wire 


6x1 


96 


Paired DB wire d 


2x4 


128 


Unpaired and paired DBs° 


4x5 


320 


Two Paired DBs d 


4x5 


320 


Paired, Unpaired and paired DBs c 


4x6 


384 


Unpaired, paired and unpaired DBs c 


4x6 


384 



a Multiples of 3.84A - lattice constant of bulk unit cell. 
^Number of Si atoms only per unit cell. 
c in [110] direction. 
d in [110] direction. 



in [110] and [110] directions as shown in Fig. 4(a) and 
(b) respectively. For unpaired DB wires, in [110] and 
[110] direction, a near-midgap state similar to that of an 
isolated unpaired DB is introduced. However, the disper- 
sion in the surface state due to an unpaired DB wire in 
[110] direction is larger than that of an isolated unpaired 
DB, whereas the dispersion in surface states due to an un- 
paired DB wire in [110] direction is only slightly larger. 
Similarly, for wires consisting of paired DBs in [110] and 
[110] direction, tt and 7r* states are introduced close to 
E v and E c respectively, which are similar to those intro- 
duced by an isolated paired DB. The dispersion in surface 
states of these wires follow similar pattern as that of un- 
paired DB wires in [110] and [110] directions respectively. 
Finally, we report electronic band structure calculations 
for isolated clusters of DBs as shown in Fig. 5, reporting 
that electronic structure of clusters can be interpreted as 
combinations of those of unpaired and paired DBs, but 
sufficiently dispersed due to interaction within the DBs. 

This paper is organized in four sections. First we talk 
about atomic structure and the assumptions made. Then 
we briefly discuss the theoretical approach. Finally, re- 
sults and conclusions are presented. 



II. ATOMIC STRUCTURE AND 
ASSUMPTIONS 

We use atomic structure of reconstructed Si(100)- 
(2x1) surface as reported by Ramstad et al 38 where the 
top four layers are relaxed due to surface reconstruction. 



FIG. 2: (color online) Ball and stick model for the unit cells of different configurations shown to demonstrate the approach 
used to calculate the electronic structure. The dimensions of unit cells (in [110] and [110] direction) for isolated DB, DB wire 
and DB cluster are adopted in such a way that the DB in one unit cell does not interfere with the neighboring unit cells in 
any inadvertent way. Note that A± and A2 are the unit vectors in [110] and [110] directions respectively, (a) Isolated unpaired 
DB on H:Si(100)-(2xl) surface, (b) Unpaired DB wire in [TlO] direction on H:Si(100)-(2x 1) surface, (c) Unpaired DB wire 
in [110] direction on H:Si(100)-(2x 1) surface, (d) Part of the surface Brillouin zone used in the calculations. Values of (fci, 
fe) at different symmetry points (r, J, K, J f ) in surface Brillouin zone are shown corresponding to [110] and [110] directions 
respectively. 



Ramstad et aJ3& report structure for five layers of Si(100)- 
(2xl)-AD surface and Si(100)-(2xl)-SD surface. We add 
eleven bulk layers to make a sixteen layer structure and 
finally passivate the bottom layer with H in dihydride 
configuration as shown in Fig. 1(a) and (b) respectively. 
We then add H atoms to each of the Si dimer atoms of 
Si(100)-(2xl)-SD structure, resulting in H:Si(100)-(2xl) 
structure as shown in Fig. 1(c). Furthermore, we remove 
one of the H atoms from H:Si(100)-(2xl) so that the re- 
sulting structure has an unpaired DB as shown in Fig. 
1(d). Fig. 1(e) shows a H:Si(100)-(2xl) surface with 
directions along and perpendicular to the dimer rows re- 
ferred to as [110] and [110] respectively. Apart from this, 
we add H atoms with H-Si bond length = 1.47 A and H-Si 
dimer angle = 111°. The bond length and the angle were 
obtained by structural optimization of a S19H14 cluster 
using local spin density approximation with 6-31 lg* ba- 
sis set as shown in Fig. 1(f) - calculations are performed 
using Gaussian 39 . The atomic co-ordinates for Si(100)- 
(2xl)-AD unit cell (Fig. 1(a)) and H:Si(100)-(2xl) unit 
cell (Fig. 1(c)) are given her o 4Q i 41 . These four unit cells, 
shown in Fig. l(a-d), are used as building blocks to con- 
struct the bigger unit cell for doing isolated DB, DB wire 



and cluster calculations. In [110] and [110] direction, the 
unit cell consists of varying layers of Si atoms depend- 
ing upon the requirement. For example, to study an 
isolated unpaired DB, we need to include enough H:Si 
layers in [110] and [110] direction so that the DBs in 
adjacent unit cells do not influence each other. We ob- 
tain this information from a previous study-i^ that a DB 
affects neighboring atoms up to approximately 10 A. Ta- 
ble I summarizes the dimensions of the unit cells and 
the number of atoms used to simulate different surfaces. 
Three representative examples of the unit cells are shown 
in Fig. 2. We show (4x4), (4x1) and (2x4) unit cells 
in Fig. 2(a-c), which are used for calculating electronic 
band structure of an isolated unpaired DB and unpaired 
DB wire in [110] and [110] directions respectively. Since, 
DBs interact electronically over a distance of approxi- 
mately 10 A, a (4x4) unit cell makes certain that DBs 
in consecutive unit cells are electronically isolated. The 
calculated electronic structure is thus that of an isolated 
DB. Similarly, DB wires and DB clusters are also spaced 
sufficiently apart. Fig. 2(d) shows part of the surface 
Brillouin zone used in the calculation (with symmetry 

points (t,j,k, j ; )aa). 
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III. THEORETICAL APPROACH 

We model the Si surface by using a periodic system 
of unit cells repeating in [110] and [110] directions. The 

lattice is defined by pAi + qA^ , where A\ and A^ are lat- 
tice unit vectors in [110] and [110] directions; p and q are 
indices representing the repeating unit cells in A\ and A2 
directions respectively. The number of neighboring unit 
cells depend on the size of the unit cell. Hamiltonian (H) 
and overlap (S) matrices are computed in EHT scheme 
using parameters developed by Cerda et al 33 . With pe- 
riodic boundary conditions in [110] and [110] directions, 
Fourier transform technique is used to transform H and S 
from the real space to the reciprocal (k) space following 
the scheme as follows, 



FIG. 3: (color online) Calculated electronic band structure 
and transmission plot for H: Si(100)-(2x 1), Si(100)-(2x 1)- 
AD and Si(100)-(2xl)-SD surface, (a) Band structure of H: 
Si(100)-(2x 1) surface shows clean band gap. All the states are 
shown for this surface. For Si(100)-(2xl)-AD and -SD sur- 
face, only 7r and 7r* states are shown. For Si(100)-(2x 1)-AD 
surface, ty state is around 0.2eV below E v and ty* state around 
0.4eV below E c as expected^. Similarly, ty and ty* states of 
Si(100)-(2x 1)-SD surface spread through the band gap. (b) 
Transmission through different surfaces in units of e 2 /h per 
spin. Complementary view of part (a) reporting transmission 
through the above-mentioned three surfaces. H:Si(100)-(2x 1) 
surface has clean band gap. Si(100)-(2x 1)-AD and -SD sur- 
face has 0.6eV and no band gap respectively due to ty and ty* 
states. The plots are shifted by 5 and 10 units along trans- 
mission axis for AD and SD surface respectively. 



N 

H(k) = H mn e a -^-^ (1) 

m—1 

N _^ _^ _^ 

S(k) = S m ne i7lM (2) 

m—l 

where k = k\ + &2, &i an d &2 are reciprocal lattice 
vectors of the surface Brillouin zone in [110] and [110] 
direction; n represents the center unit cell; m represents 
the neighboring unit cells; N is the total number of unit 
cells; d m — d n is the displacement between neighboring 
(m th ) and center (n th ) unit cell. Finally, we calculate 
energy Eigen values for different values of k along the 
surface Brillouin zone symmetry points (r, J, K, J'). 



For H:Si(100)-(2xl), Cakmak et alM- report that the 
dimer distance is 2.35A as compared to 2.23A for Si(100)- 
(2x 1)-SD as obtained by Ramstad et a& - suggesting an 
increase of approximately 5%. However, surface states af- 
ter passivation are well below E v and hence do not affect 
the band gap. It has also been reported elsewher o 11 ! 31 
that an unpaired DB introduces a structural perturba- 
tion in H:Si(100)-(2xl) structure, which is weak. We find 
that using this structure for Si dimer, the surface state is 
shifted by approximately 20meV. At room temperatue, 
this surface relaxation is anticipated to have small effect 
and hence in our calculations, we ignore it for simplic- 
ity. Furthermore, DB wires could go through Peierls dis- 
tortion. It has however been reported^ that the total 
energy gain due to such distortion is 14meV and is not 
anticipated to have a significant effect at room tempera- 
ture. Moreover, Jahn- Teller distortion in DB clusters at 
low temperature has been discussed in Ref^. The ef- 
fects due to these structural changes are also anticipated 
to be small at room temperature. Apart from this, the 
effect of dopant atoms and other defects (both surface 
and bulk) are ignored. All the atomic visualizations in 
this paper are done using GaussView^. 



IV. RESULTS 

Fig. 3(a) shows the electronic band diagram of 
H:Si(100)-(2xl) surface showing the surface and bulk 
states resulting in a clean band gap of 1.18eV. For an 
Si(100)-(2xl)-AD surface, it has been experimentally 
known that the tt state is 0.2eV below E v and tt* state 
is 0.4eV below E<£L. These tt and tt* states for Si(100)- 
(2x1)- AD surface are shown in Fig. 3(a), which are con- 
sistent with the experimental observation. The band- 
width of these states is 1.22eV and 1.26eV respectively. 
The band gap obtained in this case is 0.62eV. Similarly 
for Si(100)-(2xl)-SD surface, the tt and tt* states dis- 
perse through the whole band gap with bandwidths of 
1.39eV and 1.42eV respectively. Fig. 3(b) gives a comple- 
mentary view of Fig. 3(a) showing transmission through 
the above-mentioned three surfaces in the unit steps of 
e 2 /h. 

Fig. 4 shows calculations for an isolated unpaired DB 
and unpaired DB wire in [110] and [110] direction. An 
isolated unpaired DB introduces a near-midgap state 16 
with a very small bandwidth of 0.1 leV as shown in Fig. 
4(a). We attribute this small dispersion to the energetic 
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FIG. 4: (color online) Calculated electronic band structure of 
periodic unit cells containing isolated DBs and DB wires, (a) 
Band structure of periodic unit cells containing an unpaired 
DB and unpaired DB wires in [110] and [110] direction. The 
unpaired DB introduces a near midgap state whose dispersion 
is small. The surface states due to an unpaired DB wire in 
[TlO] direct ion have a large dispersion, due to DBs being 3.84A 
apart. For the unpaired DB wire in [110] direction, DBs being 
farther apart (7.68A), thus the states are less dispersed. All 
the bulk and surface states for structure B are shown; whereas 
for structure A and C, states inside or close to band gap are 
shown. Refer to Fig. 2 for visualization of atomic structures of 
A, B and C. (b) Electronic band structure of periodic unit cells 
containing a paired DB and paired DB wires in [110] and [110] 
direction. The paired DB introduces tt and tt* states very 
similar to the ones introduced by Si(100)-(2x 1)-AD surface, 
however the dispersion is much smaller. The DBs in paired 
DB wire in [110] are only 3.84A apart and interact in a fashion 
similar to AD surface leading to considerable dispersion of 
the DB states. Similarly, for the paired DB wire in [110] 
direction, because the DBs are now farther apart, the surface 
states appear more like those of an isolated DB pair. All the 
bulk and surface states for structure B are shown; whereas 
for structure A and C, states inside or close to band gap 
are shown. The symmetry points (r, J, K, J') along surface 
Brillouin zone for different configurations (A,B and C) are 
shown according to scale. 



isolation of this state from the conduction and the va- 
lence band states. For the DB wire in [110] direction, 
the bandwidth of the near-midgap state is 1.16eV, re- 
sulting in a very small band gap as shown in Fig. 4(a). 
We interpret this large dispersion due to the DBs be- 
ing in the close proximity of 3. 84 A and hence interacting 
strongly. For the DB wire in [110] direction, the distance 
between DBs is 7.68A. Since DBs are farther apart in 
this case, the bandwidth is 0.26eV. This bandwidth is 
much smaller than that of the wire in [110] direction and 
larger than an isolated DB as shown in Fig. 4(a). 

Fig. 4(b) presents calculation for an isolated paired 
DB and paired DB wire in [110] and [110] direction. The 
electronic structure of a paired DB is significantly differ- 
ent from an isolated unpaired DB. It introduces a bond- 



FIG. 5: (color online) Calculated electronic band structure 
of periodic unit cells containing different DB clusters, (a) 
Cluster of an unpaired and paired DBs in [110] direction, and 
of two paired DBs in [110] direction. All the bulk and surface 
states for the cluster of two paired DBs are shown; whereas 
for the cluster of unpaired and paired DBs, states inside or 
close to band gap are shown, (b) Cluster of paired-unpaired- 
paired DBs in [110] direction and of unpaired-paired-unpaired 
DBs in [110] direction. All the bulk and surface states for the 
cluster of paired-unpaired-paired DBs are shown; whereas for 
the cluster of unpaired-paired-unpaired DBs, states inside or 
close to band gap are shown. 



ing (tt) and anti-bonding (tt*) state below E v and up to 
0.28eV below E c respectively. Qualitatively, this behav- 
ior is similar to unpassivated reconstructed Si surface. 
However, dispersion of these states is much smaller. The 
bandwidth of this tt and tt* state is 0.14eV and 0.42eV re- 
spectively. The electronic band structure for a paired DB 
wire in [110] direction shows the similar tt and tt* states 
with bandwidths of 1.17eV and 1.08eV respectively. This 
large dispersion is due to strong coupling between the 
DBs in this configuration, being 3. 84 A apart. As ex- 
pected, the dispersion of tt and tt* states for DB pair 
wire in [110] direction is less than that of [110] direc- 
tion and more than that of an isolated DB pair. The 
bandwidth is 0.29eV and 0.57eV for tt and tt* state re- 
spectively. For an isolated paired DB and paired DB 
wire along [110] direction, only the state inside or close 
to band gap are shown, whereas, for the paired DB wire 
along [110] direction all the bulk and surface states are 
shown. Since the dimensions of the unit cell for config- 
urations A, B and C are different, the symmetry points 
(T, J, K, J') along the surface Brillouin zone are labeled 
accordingly. Qualitatively, the results reported are simi- 
lar to those of Watanabe et aU^ for unpaired DB wire in 
[110] direction and unpaired and paired DB wire in [110] 
and [110] directions. However, there are some differences 
in the shape of dispersion for the paired DB wire in these 
directions. 

Fig. 5 presents electronic band structure calculation 
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for DB clusters. The qualitative behavior can be at- 
tributed to a combined behavior of individual unpaired 
and paired DBs, which are 3.84A apart. Fig. 5(a) shows 
states due to a cluster of unpaired and paired DBs, 3.84A 
apart. The resulting states consist of one near-midgap 
state with bandwidth of 0.1 2eV and two tt and tt* state 
with bandwidth of 0.26eV and 0.22eV respectively. Qual- 
itatively, we interpret these states as a combination of 
states due to an unpaired and paired DB respectively. 
However, the dispersion is larger, which we attribute to 
the mutual couplings of these unpaired and paired DBs. 
Similarly, states due to a cluster of two paired DBs, 3.84A 
apart, result in two tt and two tt* states, as shown in Fig. 
5(a). These states are similar to the ones due to isolated 
paired DB. However, one of these tt states (plotted as 
green dotted line) has very small bandwidth of 83meV 
that needs to be explored further. The bandwidth of 
other tt state is 0.25eV and that of two tt* states is 0.2eV 
and 0.22eV. Fig. 5(b) extends the cluster size by incor- 
porating two clusters of (i) paired-unpaired-paired DBs, 
each 3. 84 A away, and (ii) unpaired-paired-unpaired DBs, 
each 3. 84 A away. Qualitatively, the behavior can again 
be interpreted as combination of surface states due to 
individual paired and unpaired DBs. However, we again 
find a tt state of relatively small bandwidth of 76meV 
that needs further exploration. 



V. CONCLUSIONS 

We have reported the electronic band structure calcu- 
lations showing differences between states induced in Si 



band gap due to isolated DB, DB wires and clusters. We 
extend on previous studies and report that the details 
of these states inside band gap vary considerably from 
isolated DBs to DB wires and clusters. An unpaired DB 
behaves completely differently from a paired DB. Simi- 
larly, isolated DBs are different from DB wires and clus- 
ters. Furthermore, we report that an isolated unpaired 
DB and paired DB have a very small dispersion in their 
induced states. However, wires introduce a larger dis- 
persion that can be attributed to the additional coupling 
between the DBs. Similarly, DB clusters broaden the in- 
duced states depending upon the coupling between the 
individual unpaired and paired DBs. We have ignored 
contribution of dopant atoms, dephasing processes, fur- 
ther surface relaxation effects and charge transfer effects 
in our calculations. 
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